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Abstract. This is a rather comprehensive study on the dynamics of Navier- 
Stokes and Euler equations via a combination of analysis and numerics. We 
focus upon two main aspects: (a), zero viscosity limit of the spectra of linear 
Navier-Stokes operator, (b). heteroclinics conjecture for Euler equation, its 
numerical verification, Melnikov integral, and simulation and control of chaos. 
Besides Navier-Stokes and Euler equations, we also study two models of them. 
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1. Introduction 

The dynamics of Navier-Stokes and Euler equations is a challenging problem. 
In particular, such dynamics can be chaotic or turbulent. The main challenge 
comes from the large dimensionality of the phase space where the Navier-Stokes 
and Euler equations pose extremely intricate flows. Here the dynamics we refer 
to is the so-called Eulerian, in contrast to the so-called Lagrangian, dynamics of 
fluids. The Eulerian dynamics is a dynamics in an infinite dimensional phase space 
(e.g. a Banach space) posed by the Cauchy problem of either the Navier-Stokes 
or the Euler equations as partial differential equations. The Lagrangian dynamics 
of fluid particles is a dynamics of a system of two or three ordinary differential 
equations with vector fields given by fluid velocities. The Lagrangian dynamics of 
fluid particles in three dimensions can be chaotic even when the Eulerian dynamics 
is steady (e.g. the ABC flow [2|). Nevertheless, as shown in Appendix B, the 
Lagrangian dynamics of 2D inviscid fluid particles is always integrable. 

Even though the global well-posedness of 3D Navier-Stokes and Euler equations 
is still an interesting open mathematical problem, 3D Navier-Stokes and Euler 
equations have local well-posedness which is often enough for a dynamical system 
study in the phase space. To begin such a dynamical system study, one needs to 
understand the spectra of the linear Navier-Stokes and/or Euler operators [17) . 
The spectra of the linear Navier-Stokes operators consist of eigenvalues, whereas 
the spectra of the linear Euler operators contain continuous spectra. Existence of 
invariant manifolds can be proved for Navier-Stokes equations [17) . but is still open 
for Euler equations. The size of the invariant manifolds for Navier-Stokes equations 
tends to zero in the zero viscosity limit |17j . We find that the spectra of the linear 
Navier-Stokes and Euler operators can be classified into four categories in the zero 
viscosity limit: 

(1) Persistence: These are the eigenvalues that persist and approach to the 
eigenvalues of the corresponding linear Euler operator when the viscosity 
approaches zero. 

(2) Condensation: These are the eigenvalues that approach and form a con- 
tinuous spectrum for the corresponding linear Euler operator when the 
viscosity approaches zero. 

(3) Singularity: These are the eigenvalues that approach to a set that is not in 
the spectrum of the corresponding linear Euler operator when the viscosity 
approaches zero. 
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(4) Addition: This is a subset of the spectrum of the hnear Euler operator, 
which has no overlap with the zero viscosity Hmit set of the spectrum of 
the hnear NS operator. 

We also find that as the viscosity approaches zero, the spectrum of the linear Navier- 
Stokes operator undergoes a fascinating deformation. Focusing upon the persistent 
unstable eigenvalue, we propose a heteroclinics conjecture, i.e. there should be 
a heteroclinic orbit (in fact heteroclinic cycles) associating to the instability for 
Euler equation. We will present both analytical and numerical study upon this 
heteroclinics conjecture. Then we conduct a Melnikov integral calculation along 
the numerically obtained approximate heteroclinic orbit. We also compare the 
Melnikov prediction with the numerical simulation and control of chaos for the 
Navier-Stokes equations. Numerically we mainly use the Liapunov exponent as a 
measure of chaos. In some case, we also plot the Poincare return map. We realize 
that the size of Galerkin truncations for the full Navier-Stokes equations is limited 
by the computer ability. Thus we propose two simpler models of the Navier-Stokes 
equations. For the so-called line model, we obtain numerically exact heteroclinic 
orbits for any size of Galerkin truncations. We also realize that due to viscosity, the 
chaos in Galerkin truncations of Navier-Stokes equations is often transient chaos, 
i.e. the chaos has a finite life time. Infinite life time chaos can be observed in 
Galerkin truncations of Euler equations. 

Chaos and turbulence have no good averages [21] . The matter is more funda- 
mental than just poor understanding of averages. The very mechanism of chaos 
leads to the impossibility of a good average |16j . On the other hand, chaos and 
turbulence are ubiquitous. In high dimensional systems, there exists tubular chaos 
[16j [13j |14j [15] which further confirms that there is no good average. The hope 
is that chaos and turbulence can be controlled. Two aspects of control are practi- 
cally important in applications: Taming and enhancing. When an airplane meets 
turbulence, it will be nice, safer and economic if we can tame the turbulence into 
a laminar flow or a less turbulent flow [I] [8] |25j . In a combustor, enhancing 
turbulence can get the fuel and oxidant mixed and burned more efficiently |25j [8] . 
Theoretically, one can also make use of the ergodicity of chaos to gear an orbit 
to a specific target [ 27j . Many other possibilities of applications of control can be 
designed too. An advantage of the control is that it can be done in a trial-correction 
manner without a detailed knowledge of turbulence. 

Clearly, control of chaos and turbulence has great industrial value. From a 
mathematical point of view, the question is how much mathematis is in this control 
theory. So far, the mathematical merit of the theory of control of chaos is not nearly 
as great as proving the existence of chaos [16 . Obviously, a lot of good numerics 
is in this control theory. In this article, we will address this control theory from a 
mathematical perspective, and try to formulate some good mathematical problems. 
One can add a control to any equation. But the only meaningful controls are the 
ones that are practical. Consider the 3D Navier-Stokes equations for example 

Ui^t + UjUij = -p^i + eutjj + fi + Ci , 

defined on a spatial domain T) with appropriate boundary conditions, where e = 
1/Re is the inverse of the Reynolds number, and fi = fi(t,x) is the external force. 
Assume that without the control Ci, the solutions are turbulent. The goal is to find 
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a practical control to either tame or enhance turbulence. For instance, a practical 
control Ci — Ci{t, x) should be spatially localized (perhaps near the boundary). 

Recently, there has been quite amount of works on numerical simulations of 
chaos in Navier-Stokes equations [3] [4] [7] [26] [31] [30] . Here we try to com- 
bine numerics with analysis in terms of Melnikov integrals. Unlike the sine-Gordon 
system studied in Appendix A, analytical calculation of the Melnikov integrals is 
not feasible at this moment for Navier-Stokes equations. So we will resort to nu- 
merical calculations. It is an interesting open mathematical problem that whether 
or not 2D Euler equation is integrable as a Hamiltonian system in the Liouville 
sense. Since 2D Euler equation possesses infinitely many constants of motion, it 
is tempting to conjecture that 2D Euler equation is integrable. Another support 
to such a conjecture is that both 2D and 3D Euler equations have Lax pairs [11] 
[23j [20) . In fact, it is even rational to conjecture that 3D Euler equations are 
integrable. As mentioned above, we propose the so-called heteroclinics conjecture 
for Euler equations, i.e. there exist heteroclinic cycles for Euler equations. We nu- 
merically simulate the heteroclinic orbits and use the numerical results to conduct 
numerical calculations on Melnikov integrals. In these numerical simulations, it is 
crucial to make use of known results on the spectra of linear Navier-Stokes and 
Euler operators [lOj [9] [17j . We use the numerical Melnikov integral as a tool for 
both predicting and controling chaos. As a measure of chaos, we calculate the Lia- 
punov exponents. We find that the calculated Liapunov exponents depend on the 
computational time interval and the precision of the computation. Nevertheless, as 
a measure of chaos, Liapunov exponents prove to be very robust. Since the chaos 
is often transient, we make comparison on the base of fixed time interval and fixed 
precision of computation. 

Our numerics resorts to Galerkin truncations. But Galerkin truncations are 
somewhat singular perturbations of Euler equations. Higher single Fourier modes 
have more unstable eigenvalues. Therefore, it is difficult to derive dynamical pic- 
tures for Euler equations from their Galerkin truncations. On the other hand, higher 
single Fourier modes have more dissipation under Navier-Stokes flows. So Galerkin 
truncations perform better for Navier-Stokes equations than Euler equations. To- 
day's computer ability still limits the size of the Galerkin truncations. With better 
future computer ability, Galerkin truncations will paint better and better pictures 
of Navier-Stokes and Euler equations. It seems also important to design special 
models that can picture special aspects of the dynamics of Navier-Stokes and Euler 
equations. 

2. Zero Viscosity Limit of the Spectrum of 2D Linear Navier-Stokes 

Operator 

We will study the following form of 2D Navier-Stokes (NS) equation with a 
control, 

(2.1) dtQ + n} = e[An + f{t, x) + bS{x)] , 

where i7 is the vorticity which is a real scalar- valued function of three variables t 
and X — {xi, X2), the bracket { , } is defined as 

{/,5} = (a.j)(a,,5)-(5../)(5.,<?) , 
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where ^' is the stream function given by, 

the relation between vorticity f2 and stream function 5" is, 

fl = dxiU2 - dx^ui = A* , 

and e = 1/Re is the inverse of the Reynolds number, A is the 2D Laplacian, f{t, x) 
is the external force, h5{x) is the spatially localized control, and h is the control 
parameter. We pose the periodic boundary condition 

Q.{t, xi -\-2tt,X2) ~ fl{t, xi,X2) — ^{t, xi,X2 + 27r/a), 

where a is a positive constant, i.e. the 2D NS is defined on the 2-torus T^. We 
require that / and 5 have mean zero 



Of course always has zero mean. In this case, ^E* A^^fi. 

In both 2D and 3D, the linear NS operator obtained by linearizing NS at a 
fixed point has only point spectrum consisting of eigenvalues lying in a parabolic 
region [17j . On the other hand, the corresponding linear Euler can have continuous 
spectrum besides point spectrum [17] . The spectra of the linear NS and Euler 
operators can be classified into four classes in the zero viscosity limit: 

(1) Persistence: These are the eigenvalues that persist and approach to the 
eigenvalues of the corresponding linear Euler operator when the viscosity 
approaches zero. 

(2) Condensation: These are the eigenvalues that approach and form a con- 
tinuous spectrum for the corresponding linear Euler operator when the 
viscosity approaches zero. 

(3) Singularity: These are the eigenvalues that approach to a set that is not in 
the spectrum of the corresponding linear Euler operator when the viscosity 
approaches zero. 

(4) Addition: This is a subset of the spectrum of the linear Euler operator, 
which has no overlap with the zero viscosity limit set of the spectrum of 
the linear NS operator. 

2.1. A Shear Fixed Point. For the external force / = Fcosxi [h = 0), 
= rcosxi is a shear fixed point, where T is an arbitrary real nonzero constant. 
Choose a G (0.5,0.84). There is a > such that when e > e*, the fixed point 
has no eigenvalue with positive real part, and when e € [0, e.^), the fixed point has a 
unique positive eigenvalue ^17] . Notice that this unique eigenvalue persists even for 
linear Euler (e 0). In fact, for linear Euler (e = 0), there is a pair of eigenvalues, 
and the other one is the negative of the above eigenvalue. Precise statements on 
such results are given in the theorem below. Later we will discover numerically that 
some of the rest eigenvalues of the linear Navier-Stokes operator somehow form the 
continuous spectrum of linear Euler (e = 0) as e ^ 0, while others do not converge 
to the spectrum of linear Euler (e = 0) at all [17j |18j . Using the Fourier series 
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where u^k — ^ (in fact, we always work in the subspace where aU the uj^s are real- 
valued), one gets the spectral equation of the linearized 2D Navier-Stokes operator 
at the fixed point Q = 2cosxi, 

(2.2) An^iUJn-l - + npl'^Un - An+iUn^ 



where k e ZV{0}, p = (1,0), 



A{q,r) 



1+1 — , 

^k+np' An= A{p,k + np), and 
1 











92 


r2 



(In fact, the ^n's should be counted twice due to switching q and r, but the differ- 
ence is only a simple scaling of e and A.) Thus the 2D linear NS decouples according 
to lines labeled by k. The following detailed theorem on the spectrum of the 2D 
linear NS at the fixed point = 2cosxi was proved in [17] . 

Theorem 2.1 (The Spectral Theorem [IT]). The spectra of the 2D linear NS 
operator h2.2\) have the following properties. 

(1) (Q!fc2)^ + (fci +n)2 > 1, Vn G Z/{0}. When € > 0, there is no eigenvalue of 
non-negative real part. When e = 0, the entire spectrum is the continuous 
spectrum 

—ia\k2\, ia\k2\ 

(2) k2 ~ 0, ki = 1. The spectrum consists of the eigenvalues 

A = -en^ , n E Z/{0} . 
The eigenfunctions are the Fourier modes 



+ c.c. 



ycj„p e C , ne Z/{0} 



As e — > the eigenvalues are dense on the negative half of the real axis 
{—oo, 0]. Setting e = 0, the only eigenvalue is X — of infinite multiplicity 
with the same eigenfunctions as above. 
(3) ^2 — —1; ki = 0. (a), e > 0. For any a G (0.5,0.95), there is a unique 



(2.3) 



V32 - 3a6 - 17a4 - IGa^ 



< e*(a) < 



1 



1 



1 



2(a2 + i)(q,2 +4) *^ ^ (q,2 

where the term under the square root on the left is positive for a S 
(0.5,0.95), and the left term is always less than the right term. When 
e > e*(a), there is no eigenvalue of non-negative real part. When e — 
e*{a), X = is an eigenvalue, and all the rest eigenvalues have negative 
real parts. When e < e*(a), there is a unique positive eigenvalue A(e) > 0, 
and all the rest eigenvalues have negative real parts. e~^A(e) is a strictly 
monotonically decreasing function of e. When a G (0.5,0.8469), we have 
the estimate 



< 



2{o? + 1) ^ 4(a2 + l)(a2 + 4) 
\o?[\~o?) 2 



- t(o? + 1) < A(e) 



2(a2 
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where the term under the square root on the left is positive for a G 
(0.5,0.8469). 

Y 2(«2 + 1) 4(^2 + i)(^2 + 4) =^ ^(^) =^ y 2(^2 + 1) ■ 

In particular, as e 0^, A(e) = C(l). 

('fej. e = 0. When a € (0.5,0.8469), iwe /laue only two eigenvalues Xq 
and — Aq, where Xq is positive, 

/a2(l-Q2) q4(c,2 _^3) /a2(l-a2) 

Y 2(a2 + l) " 4(a2 + i)(^2+4) < ^ < y 2(^2 + i) • 

The rest of the spectrum is a continuous spectrum [—ia, ia] . 
(c). For any fixed a e (0.5, 0.8469), 

(2.4) lim A(e) = Aq . 

(4) Finally, when e = 0, the union of all the above pieces of continuous spectra 
is the imaginary axis iM. 

Remark 2.2. In the current periodic boundary condition case, viscosity does 
not destablize the flow in contrast to the non-sUp boundary condition case [24j . The 
Orr-Sommerfeld equation and Rayleigh equation have the same periodic boundary 
condition in the former case, and different number of boundary conditions in the 
latter case. 

The following invariant manifold theorem of the 2D NS at the fixed point 
= 2cosxi was also proved in [17] . 

Theorem 2.3 (Invariant Manifold Theorem [l7j). For any a G (0.5, 0.95), and 
e e (0, £*(«)) where e»(a) > satisfies i2.3]} . in a neighborhood of fl = 2cosa;i in 
the Sobolev space H^{T'^) (£ > 3), there are an 1-dimensional C°° unstable manifold 
and an l-codimensional C°° stable manifold. 

One of the goals of the work |17| is to study the zero viscosity limit of the 
invariant manifolds of the 2D NS. For this study, it is crucial to understand the 
deformation of the linear spectra as e ^ 0+. Below we will study this numerically. 

When fci = and k2 = 1, a — Q.7, the unique e^, in (|2.3p belongs to the interval 
0.332 < e* < 0.339, such that when e < e*, a positive eigenvalue appears. We test 
this criterion numerically and find that it is very sharp even when the truncation of 
the linear system (|2.2|) is as low as |n| < 100. As e — > 0^", we tested the truncation 
of the linear system (|2.2p up to |n| < 1024 for a = 0.7, the patterns are all the 
same. Below we present the case \n\ < 200 for which the pattern is more clear. 

Figure [T] shows the case e = 0.14 where there is one positive eigenvalue and all the 
rest eigenvalues are negative. Figure [2] shows the case e — 0.13 where a pair of 
eigenvalues jumps off the real axis and becomes a complex conjugate pair. Figure 
[3] shows the case e = 0.07 where another pair of eigenvalues jumps off the real axis 
and becomes a complex conjugate pair. Figure |4] shows the case e = 0.03 where 
another pair of eigenvalues jumps off the real axis and becomes a complex conjugate 
pair, while the former two pairs getting closer to each other. Figure [5] shows the 
case e = 0.0004 where many pairs of eigenvalues have jumped off the real axis and 
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Figure 1 . The eigenvalues of the hnear system (12.21) when fci = 
and k2 = 1, a = 0.7, and e = 0.14. 

a bubble is formed. Figure [6] shows the case e = 0.00013 where the bubble has 
expanded. As e — s- the limiting picture is shown in Figure [71 Setting e = 0, 
the spectrum of the line ki = and ^2 = 1 of the linear Euler operator is shown 
in Figure [H where the segment on the imaginary axis is the continuous spectrum. 
Comparing Figures [7] and [SI we see that the two eigenvalues represent "persistence" , 
the vertical segment represents "condensation", and the two horizontal segments 
represent "singularity". Next we study one more line: fci = and fc2 = 2 (a = 0.7). 
In this case, there is no unstable eigenvalue. Figure [9| shows the case e = 1.5 where 
all the eigenvalues are negative. As e is decreased, the eigenvalues go through the 
same process of jumping off the real axis and developing a bubble. Figure [TOl shows 
the case e = 0.00025 where the bubble has expanded. As e 0+, the limiting 
picture is similar to Figure [7| except that there is no persistent eigenvalue. The 
cases fci = and ^2 > 2 (a = 0.7) are all the same with the case fci = and fc2 = 2 
(a = 0.7). Figure [TT] shows the limiting picture of the entire spectrum of the linear 
NS operator as e ^ 0+ . Figure [T^j shows the entire spectrum of the linear Euler 
operator (e = 0). 

The fascinating deformation of the spectra as e ^ 0+ and the limiting spectral 
picture clearly depict the nature of singular limit of the spectra as e ^ O'*' . In the 
"singularity" part of the limit, there is a discrete set of values for the imaginary 
parts of the eigenvalues, which represent decaying oscillations with a discrete set of 
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Figure 2 . The eigenvalues of the hnear system (12.21) when fci = 
and ^2 = 1, a = 0.7, and e = 0.13. 



frequencies. Overall, the "singularity" part represents the temporally irreversible 
nature of the e — s- 0^ limit, in contrast to the reversible nature of the linear Euler 
equation (e = 0). 

2.2. A Cat's Eye Fixed Point. In this subsection, we will study another 
important fixed point - a cat's eye fixed point. The periodic domain now is the 
square, i.e. a = I (instead of 0.7). The cat's eye fixed point in physical variable is 
given by 

(2.5) r2 = 2cosmxi + 27C0smx2 , 

where m is a positive integer, and 7 e (0,1]. In terms of Fourier modes: Let 
p = (m, 0) and q = (0, m), then the cat's eye is given by 

and all other w^'s are zero. The spectral equation for the linear 2D NS operator at 
the Cat's eye is then given by 

Awfe = A(j),k-p)ujk-p'' A{p,k+p)u}k+p~ e\k\'^ujk 

(2.6) +^A{q, k - q)ujk-q - lA{q, k + q)ujk+q , 
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Figure 3. The eigenvalues of the hnear system (12.21) when fci = 
and k2 = 1, a = 0.7, and e = 0.07. 



where 



A{k,r) = 











k2 


f2 



First we study the case of m = 1 and 7 = 0.5. Changing the value of 7 does 
not affect the deformation patterns of the eigenvalues of the linear NS (|2.6p as 
e ^ 0"*". We truncate the system (I2.6P via the Galerkin truncation |fci| < 32 and 
1^2! < 32. This is the largest Galerkin truncation that we are able to compute in 
a reasonable time. For smaller Galerkin truncations, the deformation patterns of 
the eigenvalues are the same. When e = 150, all the eigenvalues of the linear NS 
(|2.6p are negative as shown in Figure [121 When e is decreased to e = 120, a pair 
of eigenvalues jumps off the real axis as shown in Figure [M) When e is decreased 
to e = 80, three pairs of eigenvalues jump off the real axis as shown in Figure [TSl 
When e is decreased to e = 10, many pairs of eigenvalues have jumped off the real 
axis and form several parabolas as shown in Figure [HI When e is decreased to 
e = 1 , many parabolas are formed as shown in Figure 1171 The e — > O"*" limiting 
picture of the eigenvalues of the linear NS (|2.6p is that the eigenvalues are dense 
on the entire left half plane as shown in Figure [TH] The continuous spectrum of 
the linear Euler, i.e. e = in p.6|) . in any Sobolev space iJ''(T^) where s is a 
non-negative integer, is a vertical band of width 2sa symmetric with respect to 
the imaginary axis {A : |Re(A)| < sa} as shown in Figure [HI where cr > is the 



NAVIER-STOKES AND EULER EQUATIONS 



11 



0.4r 

0.3- 
0.2- 



0.1 - 




-0.1 - 
-0.2- 
-0.3- 



-1400 -1200 -1000 -800 -600 -400 -200 200 

Re(l) 

Figure 4. The eigenvalues of the hnear system (12.21) when fci = 
and k2 = 1, a = 0.7, and e = 0.03. 

largest Liapunov exponent of the vector field given by the cat's eye (|2.5I) ,29 . Thus 
the width of the vertical band is proportional to the scale s of the Sobolev space 
i?'*(T^). The union of all such bands for all integers s > is the entire complex 
plane. The eigenfunctions of the linear NS (|2.6p when e > belong to H'^{T'^) 
for all integers s > 0. All the eigenvalues of the linear NS (|2.6p condense into the 
entire left half plane - "condensation" . The right half plane (or right half of the 
vertical band corresponding to 7J^(T^)) represents "addition". Thus the possible 
instability hinted by the right half band of the continuous spectrum of linear Euler 
in H^{T'^) can not be realized by real viscous fluids. 

Next we study the case of m = 2 and 7 — 0.5. Changing the value of 7 does not 
affect the deformation patterns of the eigenvalues of the linear NS (|2.6[) as e 0+. 
We truncate the system (|2.6p via the Galerkin truncation \ki \ < 32 and \k2\ < 32. 
We increase the value of e up to 2 x 10^, there are still eigenvalues with nonzero 
imaginary parts. These eigenvalues seem always complex no matter how large is 
e. The imaginary parts of these eigenvalues are unchanged between e = 2 x 10^ 
and e = 800 as can be seen from Figures [201 and [^T] Decreasing e, the deformation 
patterns are similar to those of m = 1. When e = 0.2, many eigenvalues have 
jumped off the real axis and form a dense parabolic region as shown in Figure [22l 
Decreasing e further, 6 eigenvalues with positive real parts appear, two of which 
are real, and the rest four are complex. The limiting picture of the spectrum of 
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Figure 5 . The eigenvalues of the hnear system (12.21) when fci = 
and fca 1, a = 0.7, and e ^ 0.0004. 



the hnear NS operator (|2.6p as e ^ is the same with the m = 1 case as shown 
in Figure [THl except the extra six unstable eigenvalues. The continuous spectrum 
of the linear Euler operator (|2.6p where e = is the same with the m — 1 case as 
shown in Figure 1191 There is no analytical result on the eigenvalues of the linear 
Euler operator (|2.6p where e = 0. The numerics indicates that the 6 eigenvalues in 
the right half plane and other 6 eigenvalues in the left half plane, of the linear NS 
operator (|2.6p persist as e ^ 0, and result in 6 eigenvalues in the right half plane 
and their negatives for the linear Euler operator (|2.6p where e = 0. 



3. The Heteroclinics Conjecture for 2D Euler Equation 

Setting e = in the 2D Navier-Stokes equation ()2.ip . one gets the corresponding 
2D Euler equation for which one has the following constants of motion: 

/ \u\^dx , / F{n)dx 

where F is an arbitrary function. Consider the simple fixed point f2 = Fcosxi 
(r ^ real constant). It has one unstable and one stable real eigenvalues which 
are negative of each other. The rest of the spectrum is the entire imaginary axis 
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Figure 6. The eigenvalues of the hnear system (12.21) when fci = 
and fca = 1, a = 0.7, and e = 0.00013. 



which is a continuous spectrum [lOj |17j . We wiU use the constant of motion 

[ n^dx- [ \u\'^dx 

to build a Melnikov integral for the corresponding 2D Navier-Stokes equation (|2.ip . 
We will try to make use of the Melnikov integral as a measure of chaos and to 
conduct a control of chaos, around the line of fixed points 57 = F cos parametrized 
by F. G is a linear combination of the kinetic energy and the enstrophy. The 
gradient of G in is given by 

which is zero along the line of fixed points H. = Fcosa;i. We define the Melnikov 
integral for the 2D NS ((2T|) as 

(3.1) M = r f VnG[An + f{t,x) + bS{x)]dxdt 

= Mo + bMc , 
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Figure 7. The (e 0+) limiting picture of the eigenvalues of the 
linear system (|2.2|) when fci = and k2 = 1, and a = 0.7. 



where 

Mo = ^ / / (n + A-^n)\An + fit, x)]dxdt , 
47r^ Jt-2 

a f 

Mc = — TT {^ + A-^n)d(x)dxdt . 

The question is: Where do we evaluate Ml We propose the following conjecture. 
• The Heteroclinics Conjecture: There is a heteroclinic orbit of the 2D Eulcr 
equation that connects f2 = Fcosxi and — f2. 
The rationality of this conjecture has been discussed in the Introduction. If this 
conjecture is true, we can evaluate M along the heteroclinic orbit. Also, under 
the perturbation of the e term, the heteroclinic orbits may break and re-connect 
somewhere, thereby generating the heteroclinic chaos. In Appendix A, we show 
that one can use a Melnikov integral as a criterion to rigorously prove the existence 
chaos, and to conduct control of chaos. In the current case of 2D NS, the phase 
space is much more complicated. We do not expect the above Melnikov integral to 
be a rigorous measure of chaos, rather we believe that it is still relevant to chaos 
and control of chaos. In fact, often the chaos in 2D NS is a transient chaos. In 
general, the zeros of a Melnikov integral do not immediately imply the existence of 
heteroclinic or homoclinic chaos, rather they imply the intersections of the broken 



NAVIER-STOKES AND EULER EQUATIONS 



15 



Figure 8. The spectrum of the Unear system (j2.2[) when e = 0, 
ki = and k2 = 1, and a = 0.7. 

heterochnic orbit with certain center-stable manifold [16. . In fact, the Melnikov 
integral represents the leading order distance between the broken heterochnic orbit 
and the center-stable manifold [16] . Here for 2D NS, rigorous justification even 
on the above claim is an open problem due to the fact that existence of invariant 
manifolds for 2D Euler is an open problem. 

The 2D Euler equation has several symmetries: 

(1) fl{t, Xl,X2) > n{t,-Xl,-X2), 

(2) Cl{t, Xi,X2) > ~Q{~t,Xi,X2), 

(3) ^{t,Xi,X2) > -fl{t, -Xi,X2), or Q{t,Xi,X2) — > -n{t,xi,-X2), 

(4) n{t,xi,x2)-^n{t,xi + ei,x2 + 02), yoi,02. 

The first symmetry allows us to work in an invariant subspace in which all the lu^ 's 
are real- valued. This corresponds to the cosine transform in l|3.2p . In this paper, we 
will always work in the invariant subspace where all the w^'s are real- valued. The 
second symmetry maps the unstable manifold of T cos xi into the stable manifold 
of — rcosxi. The third symmetry maps the unstable manifold of rcosxi into the 
unstable manifold of — rcosxi. By choosing 9i = tt, the fourth symmetry maps the 
unstable manifold of Fcosxi into the unstable manifold of — Fcosxi. To maintain 
the cosine transform, the 0i and 02 in the fourth symmetry can only be tt and n/a. 

If the heteroclinics conjecture is true, there is in fact a pair of heterochnic cycles 
due to the above symmetries. Indeed, if there is a heterochnic orbit asymptotic to 
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Figure 9. The eigenvalues of the hnear system (12.21) when fci = 
and k2 ~ 2, a = 0.7, and e ~ 1.5. 
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Fcosxi and — Fcosxi as i — > — oo and +cx), then the third symmetry generates 
another heteroclinic orbit asymptotic to — F cos xi and F cos xi as i — oo and +oo. 
Together they form a heteroclinic cycle. Finally the second symmetry generates 
another heteroclinic cycle. That is, we have a pair of heteroclinic cycles. 
Using the Fourier series 

(3.2) J2 Wfce*('=i^i+"'=^^^) , 

fcGZV{0} 

where uj^k = ^ and F-k = Fk, one gets the kinetic form of the 2D Euler equation 

UJk= ^ A{m, n) UJrnl^n , 
k—m-\-n 

where 

mi ni 

7712 «2 

Denote by E the hyperplane 

T, — {lo \ LUk = , V even fc2} • 

Notice that the existence of invariant manifolds around the fixed point ~ Fcosxi 
is an open problem. We have the following theorem. 



A{m, 77 ) 



77? + (Q!772)^ TR? + (q;7772)^ 
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Figure 10. The eigenvalues of the hnear system (|2.2p when fci = 
and fca = 2, a = 0.7, and e ^ 0.00025. 



Theorem 3.1. Assume that the fixed point D, = Tcosa;! has a 1- dimensional 
local unstable manifold , and VF" HE 7^ 0; then the heteroclinics conjecture 
is true, i.e. there is a heteroclinic orbit to the 2D Euler equation that connects 
51 — rcosa;i and — f2. 

Proof. Let il{t, xi,X2) be an orbit in parametrized such that 
Then by the definition of E, 

(3.3) fl{0,Xi, X2) = — ri(0, Xi, 0:2 + 7r/a) . 

By the second and fourth symmetries, 

-ri(-i, xi,X2 + vr/a) 

is in the stable manifold of —51. Thus 

f2(t, xi, 2:2) and — fl{—t, xi, X2 + tt/q:) 

are connected at t = 0, and together they form a heteroclinic orbit that connects 
= Fcosa;! and —il. □ 
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Figure 11. The (e 0+) limiting picture of the entire spectrum 
of the linear NS operator (|2.2|) when a ~ 0.7. 



4. Numerical Verification of the Heteroclinics Conjecture for 2D Euler 

Equation 

Besides the symmetries mentioned in last section, we will also make use of the 
conserved quantities: kinetic energy E — ^ \k\~'^Lu'^ (where j/cp — kf + a^fc^) and 
enstrophy S — '^uj'^, which will survive as conserved quantities for any symmetric 
Galerkin truncation, to help us to track the heteroclinic orbit. We will only consider 
the case that all the w^'s are real- valued (i.e. cos-transform). 

We make a Galerkin truncation by keeping modes: {|fci| < 2, \k2\ < 2}, which 
results in a 12 dimensional system. We choose a — 0.7. After careful consideration 
of the above mentioned symmetries and conserved quantities {E = S = 1), we 
discover the following initial condition that best tracks the heteroclinic orbit: 

^U,o) = ^u,2) = , V j , 
(4.1) W(oa) = 0.603624 , = -w(-i,i) = 0.357832 , 

a;(2,i) = ^{-2,1) = 0.435632 . 

Starting from this initial condition, we calculate the solution in both forward and 
backward time for the same duration of T = 11.8, and we discover the approximate 
heteroclinic orbit asymptotic to 2cosa;i and — 2cosa;i as t —>■ —oo and -l-oo, as 
shown in Figure [23l Then the third symmetry generates another heteroclinic orbit 



NAVIER-STOKES AND EULER EQUATIONS 



19 



Figure 12. The spectrum of the hnear Euler operator 
6 = and a = 0.7. 



where 



asymptotic to — 2cosa::i and 2cosa;i as i — > — oo and +oo. Together they form 
a heterochnic cycle. FinaUy the second symmetry generates another heterochnic 
cycle. That is, we have a pair of heterochnic cycles. Notice also that the approxi- 
mate heterochnic orbit in Figure[23]has an extra loop before landing near —2 cosa;i. 
This is due to the k2 — 2 modes in the Galerkin truncation. For smaller Galerkin 
truncations, the heterochnic orbits can be calculated exactly by hand and have no 
such extra loop |12| [22] , and existence of chaos generated by the heterochnic orbit 
can be rigorously proved in some case 



Remark 4.1. We have also conducted numerical experiments on Galerkin trun- 
cations by keeping more modes: {|fci| < 4, \k2\ < 4} and {|fci| < 8, |fc2| < 8}. We 
found orbits that have similar behavior as the approximate heterochnic orbit in 
Figure [221 but their approximations to hetcroclinics are not as good as the one in 
Figure m 

5. Melnikov Integral and Numerical Simulation of Chaos in 2D 
Navier-Stokes Equation 

Without the control (6 = 0), using Fourier scries for the 2D NS equation (12.11) . 
feezVfo} fcezVfo} 
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Figure 13. The spectrum of the hnear NS operator 
e = 150, m — I and 7 = 0.5. 



where 



where W-fe = LUk and F-k = Fk (in fact, we always work in the subspace where all 
the Wfc's and Ffc's are real- valued) , one gets the kinetic form of 2D NS 

Wfc = ^ A{m,n) LLijni^ri + f {-[kl + {ak2)'^]uJk+Fk) , 

k—m-\-n 



where 



A{m, n) 



mi ni 
m2 rL2 



For the numerical simulation of chaos, we continue the study on the Galerkin 
truncation: {|fci| < 2, \k2\ < 2}. We will use the Melnikov integral ()3.1|) to test the 
existence of chaos. We always start from the initial condition (|4.ip . We choose the 
external force 

(5.1) / = asini cos(a;i + 01x2) . 
Then the Melnikov integral (13. ip has the expression 

(5.2) Mo = Ml + a y^M| + M| sin(to + 6*) , 
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Figure 14. The spectrum of the hnear NS operator 
e = 120, TO = 1 and 7 = 0.5. 



where 



where 

Ml 
Ah 
M3 



M: 



cost 



/+00 /'27t/o /'27T 
/ {^ + A-^n)An dxidx2dt , 

^00 Jo Jo 

/-\-OD i>2-K / a i>2-K 
/ A^-^il) cosi cos(xi + ax2) dxidx2dt , 

- 00 1/ '-'0 

— TT / / / (Q + A~^fl)smtcos(xi + 0x2) dxidx2dt , 
47r^ J-00 Jo Jo 



where fl{t) is the approximate heterochnic orbit in Figure [23] with Q{0) given 
by (|4.ip . The time integral is in fact over the interval [—11.8,11.8] rather than 
(—00,00). Direct numerical computation gives that 

Ml = -29.0977 , M2 = -0.06754695 , 1/3 = 0. 

Setting Mo = in (|5.2p . we obtain that 

. , , 430.77741 
sm(io + tt) = . 
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Figure 15. The spectrum of the hnear NS operator 
e = 80, m = 1 and 7 0.5. 



where 



Thus, when 

(5.3) \a\ > 430.77741 , 

there are solutions to Mq = 0. Next we will test the Melnikov criterion (|5.3p and see 
if it is related to chaos. We define an average Liapunov exponent a in the following 
manner: For a large time interval t g [0, T], let to = T and 

tn^T + n2TT , where < n < iV , and TV = 10^ or 2 x 10^ . 

We define 



A In 
27r 



\ALu{tn + 2n)\\ 
||Aa;(i„)ll 

Then the average Liapunov exponent a is given by 

^ = ]^ E ^" ■ 

n=0 

We introduce the Poincare return map on the section given by Ci;(l,0) = and 
we only record one direction intersection (from a;(l,0) positive to negative). For a 
large time interval t € [0,T], we only record the last 1000 intersections and we use 
• to denote the intersections with the two heteroclinic cycles. All the numerical 
simulations start from the initial condition (|4.ip . The average Liapunov exponent 
computed here depends on the time interval, the ensemble of average, and computer 
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Figure 16. The spectrum of the hnear NS operator (|2.6p where 
e = 10, m = 1 and 7 = 0.5. 



accuracy. It only makes sense when it is compared in the same setting. When e = 0, 
there is no dissipation and no forcing. For a large time interval t S [0, T], the average 
Liapunov exponent a is as follows: 

T = 4 X IO^'tt T = 8 X lO^TT T = 12 x IO^'tt T = 80 x IO^'tt , 
a = 0.042 a = 0.0344 a = 0.044 a = 0.0848 . 

Figures 12411251 are the corresponding Poincare return map plots. The dynamics is 
chaotic. It seems that the life time of the chaos is infinite (i.e. non-transient chaos). 

When e > 0, we find that in all the cases that we tested, the dynamics is 
always a transient chaos. The Melnikov criterion is only some sort of necessary 
condition for the existence of heteroclinic chaos |16j . When the Melnikov integral 
is zero, it gives an indication of a re-intersection of the broken heteroclinic orbit 
with certain large dimensional center-stable manifold [16j . We believe that such 
a re-intersection will be reflected by the Liapunov exponent as inducing transient 
chaos. When e = 10-^ a e [0, 1208], and T = 4 x IO-^tt, we find that 
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Figure 17. The spectrum of the hnear NS operator 
e ~ 1, m ~ 1 and 7 ~ 0.5. 



where 



For instances, 
(5.4) 

a = 400 a = 430 a = 440 a = 650 a = 850 , 

cr = 4.9x10-4 cr = 2.6x10-* (7 = 5.3x10-4 a = 5.9 x 10"^ cr = 1.0 x lO"* 

But we discover a sharp jump of a near a = 1208 as shown below: 
(5.5) 

T = 2 x lOV T = 4 X lO^TT T = 8 X IO^tt T = 8 x IO'^tt , 
a = 1208 (7 = 3.6 X lO""* ct = 3.9 x 10""* ct = 4.1 x lO""* cr = , 
a = 1208.2 cr = 6.1x10^2 cr = 6.1 x 10"^ cr = 2.5 x 10"^ cr = . 

For T = 4x lO^TT, the corresponding Poincare return map plots are shown in Figures 
[2611271 When T = 80 x IO^'tt, the chaos for the a = 1208.2 case also disappears. 
When a > 1208.2, a can still be ~ lO"''. But we did not observe any sharp jump of 
(J. We believe that the sharp jump of a near a ~ 1208 is due to the generation of 
new transient heteroclinic chaos. Our Melnikov integral calculation predicts when 
\a\ > 430.77741, the broken heteroclinic orbit re-intersects with certain center- 
stable manifold. This does not mean that new heteroclinic orbit is automatically 
generated. Perhaps near a = 1208, new heteroclinic cycle is generated and leads to 
transient heteroclinic chaos. Since our analysis cannot access such information in 
the phase space, our comments here are purely speculations. 
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Figure 18. The limiting picture of the spectrum of the hnear NS 
operator (|2.6p as e ^ 0, where m = 1 and 7 = 0.5. 



Remark 5.1. We have also conducted numerical experiments on Galerkin trun- 
cations by keeping more modes: {|fci| < 4, |fc2| < 4} and {\ki\ < 8, |fc2| < 8}. We 
found that when e = 0, the strength of chaos increases as the modes are increased: 
For T = 8 X lO^TT, 

|fcl|,|fc2|<2 |fci|,|fc2|<4 |fci|,|fc2|<8, 

cr = 4.7x10^2 (7 = 1.3x10^1 (7 = 1.7x10^1. 

Nevertheless, this does not hint that the dynamics of 2D Euler equation is chaotic 
since all Galerkin truncations are perturbations of the 2D Euler equation. In fact, 
higher single Fourier modes (as fixed points) have more eigenvalues with positive 
real parts. 

When e > 0, the strength of chaos decreases as the modes are increased. Higher 
modes have more dissipations. Also all the chaos are transient. After enough time 
(~ 2 X IO^tt), the e = chaos is almost smeared away by dissipation, we believe 
that at this stage the re-intersected heteroclinic orbits play a role and can enhance 
the transient chaos. It is this stage where the Melnikov calculation is effective. 




6. Melnikov Integral and Control of Chaos in 2D Navier-Stokes 

Equation 

Now we turn on the control (6 7^ 0). we eontinue the study on the Galerkin 
truncation: {|fci| < 2, |fc2| < 2}. We choose S{x) as foUows 



(6.1) ^(x) ^^e^^^^i^i 



+ak2X2) 



Then the Melnikov integral M in (13. ip is given by 
(6.2) M = Ma + bMc , 

where Mq is given by (|5.2p and 



/+00 p27T / a n27T 
/ {^ + A-^n)6{x) dxidx2dt 

-00 Jq Jo 



evaluated along the approximate heteroclinic orbit in Figure [23l We find that 

Mc = -18.6884 . 

Thus 

M = -29.0977 - 18.6884 b + 0.06754695 a sin(io + tt) . 
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Figure 20. The spectrum of the hnear NS operator 
e = 20000, m = 2 and 7 = 0.5. 



where 



When 

6:= -1.557 , 

the Melnikov integral M has roots for any a ^ Q. 

All the numerical simulations start from the initial condition (|4.ip . When 
e = 10-5^ b = -1.557, and T = IO^tt, we find that: 



a = 


1 


a = 


10 


a ~ 


200 


a ~ 


400 


a 


= 800 , 


a = 


7.1 X IQ-^ 




8.8 X 10^4 


a = 


9.3 X 10-4 


a = 


9.0 X 10-4 


(7 


6.6 X 10-4 


a = 


1000 


a = 


1208 


a = 


1208.2 


a = 


1500 


a 


= 3000 , 


a = 


8.6 X 10--* 


a = 


8.1 X 10-* 


a = 


8.4 X IQ-* 


a ~ 


8.6 X 10-4 


a 


= 7.7 X 10-4 



In comparison with (|5.4p . the values of the Liapunov exponents under the control 
are doubled. Thus the control seems enhancing chaos but not dramatically. We did 
not observe the sharp jump of the values of a around a = 1208 as in the 5 = case 
(|5.5p . Here the control theory is not as rigorous as that of sine-Gordon system in 
Appendix A. But we believe that the Melnikov integral can play a significant role 
in the control of chaos in NS equation. After all, the chaos in NS is generated by 
instabilities characterized by unstable eigenvalues. And these unstable eigenvalues 
persist for Euler equation as shown in a previous section. For Euler equation, 
these unstable eigenvalues characterize hyperbolic structures which are very likely 
of heteroclinics type due to infinitely many constants of motion. Thus Melnikov 
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Figure 21. The spectrum of the hnear NS operator (|2.6p where 
e = 800, TO = 2 and 7 = 0.5. 



integrals supported upon these hyperbolic structures should play an important role 
in predicting and controling chaos. 



7. Zero Viscosity Limit of the Spectrum of 3D Linear Navier- Stokes 

Operator 

We will study the following form of 3D Navier-Stokes equation with a control, 

(7.1) dtQ + {u ■ \7)Q - {Q ■ V)u = e[AO + f{t, x) + bS{x)] , 

where u = (mi,U2;U3) is the velocity, Q — (fii, 172, rJa) is the vorticity, V — 
{dxi,dx2,dx3), fl — Wxu, W-u~Q,e — 1/Re is the inverse of the Reynolds 
number, A is the 3D Laplacian, and f{t,x) = (fi(t, x), f2(t, x), f3{t, x)) is the ex- 
ternal force, bS{x) is the spatially localized control, and b is the control parameter. 
We also pose periodic boundary condition of period (27r/a, 27r//3, 27r), i.e. the 3D 
NS is defined on the 3-torus T^. We require that u, fl, / and 6 all have mean zero. 
In this case, u can be uniquely determined from Q by Fourier transform: 

Ui{k) = l\k\-^[k2LU3{k) ~ ksiO^ik)] , 

U2{k) = i\k\-^[k3u;i{k) ~ kitusik)] , 
Usik) = i\k\-^[kiu;2{k) - k2cui{k)] , 
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Figure 22. The spectrum of the hnear NS operator (|2.6p where 
e = 0.2, TO = 2 and 7 = 0.5. 



which can be rewritten in the compact form 

Uiik) = i|fc|"^Q„„fc„a;„(fc) , 
where Simn is the permutation symbol {i, to, n — 1, 2, 3), 

k = {ki,k2,h) = {aKi,l3K2,K3) , « = (k1,K2,K3) , 

Using these Fourier transforms together with 

ksz3/{o} Kez3/{o} 
we can rewrite the 3D NS (17.11) into the kinetic form 



ik-x 



dtuji{k) + ks ^ ^kmUJn{k)[e£„inUJs{k) - esmn(^i{k)] 



k=k + k 



(7.2) 



\k\^iOe{k) + FAk) + bAe{k)] 



A popular example of fixed points of the 3D NS (|7.ip is the so-called ABC flow [2] 
(7.3) ui — A sin 2:3 + C cos X2, U2 = B sin xi + A cos 2:3 , U3 = C sin a;2 + i? cos xi , 
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Figure 23. The approximate heteroclinic orbit projected onto the 
(a;(i_o), ^(i,i))-plane in the case of the {|fci| < 2, \k2\ < 2} Galerkin 
truncation of the 2D Euler equation. 



where a = (3 — 1 and Q = u = f and 6 = 0. The popularity comes from the fact 
that the Lagrangian fluid particle flow generated by the vector field (|7.3p can still 
be chaotic [2]. On the other hand, in Appendix B, we prove that the Lagrangian 
flow generated by any solution to the 2D Euler equation is always integrable. 

7.1. A 3D Shear Fixed Point. Below we will study the simplest fixed point 
- the 3D shear flow (which is also a special case of the ABC flow (|7.3p where A = 2 
and B = C = 0): 

(7.4) fii — 2 sin 2:3 , ^2 — 2cosx3 , = . 

Let p — (0, 0, 1), the Fourier transform lu* of the flxed point is given by: 

and uig{k) = 0, Vfc 7^ p or —p. We choose a = 0.7 and /3 = 1.3 hoping that the 
flxed point uj* has only one unstable eigenvalue. The spectral equations of the 3D 
linear NS operator at the flxed point w* are given by 

[{ki + ik2) - ik2\k - p\^'^]ijJi{k ~ p) 

+ [i\k~p\-'^ki+i\k-p\-^{ki +ik2){kj, - \)]uj2{k-p) 

+ - ik2{ki + ik2)\k - p\-^]uj3{k -p) + €\k\'^uoi{k) 

+ [-{ki - ik2) - ik2\k + p\^'^]uji{k + p) 

+ [i\k + p\-'^ki - i\k + p\-^{ki ~ ik2)ik3 + l)]uj2{k+p) 

+ [-1 + ik2iki ~ ik2)\k + p\-'^]uj3{k + p) = -Awi(fc) , 
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Figure 24. The Poincare return map plot projected onto the 
plane (cjjo.i), Ci;(i^i))-plane, in the case of the < 2, |fc2| < 2} 

Galerkin truncation of the 2D Euler equation (i.e. e = 0), where 
t € [0,T], T = 4 X IO^'tt, only the last 1000 intersections are 
recorded. 



[fcalfc - p\-^ - i{h - l)(fci + ik2)\k - pr']wi(fc - p) 

+ [{ki + ik2) - ki\k - p\-^]L02{k - p) 

+ + iki{ki + ik2)\k - p\-^]uo:i{k - p) + e|fc| 

+ [-fc2|fc +pr^ + i(fc3 + l)(fci - ik2)\k + p\-'^]uji{k + p) 

+ [-{ki - ik2) + ki\k + p\-'^]uj2{k + p) 

+ [i - iki{ki - ik2)\k + p\^'^]u;3{k + p) = -\uj2{k) , 

[ik2{ki + ik2)\k - p) + [~iki{ki + ik2)\k ~ p\^'^]uj2{k - p) 

+ [ki + ik2]uj3{k — p) + e|fcpCLi3(fc) + [—ik2{ki — ik2)\k + p|^^]cji(fc + p) 
+ [iki{ki - ik2)\k + p\^^]uj2{k + p) + [-{ki - ik2)]uJ3{k + p) = -Xujsik) ■ 

Thus the 3D linear NS operator also decouples according to the lines k+jp {j € Z). 
Next we study the zero viscosity limit of the spectrum of this 3D linear NS operator. 
When k = (a, 0, 0), we have tested the truncation of the hue k + jp up to \j\ < 400. 
The deformation pattern stays the same. Below we present the case |j| < 100 for 
which the pictures are more clear. Figure [25] shows the case e = 2.0 where all 
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Figure 25. The same as Figure [Ml except T = 80 x IO^tt. 



the eigenvalues are negative. Figure [29] shows the case e = 1.8 where a pair of 
eigenvalues jumps off the real axis. When e < 0.66, a unique positive eigenvalue 
appears. Figure [281 shows the case e = 0.0007 where a bubble has developed. 
As e ^ 0+, the limiting picture is the same with Figure [71 When e = 0, The 
spectrum picture is the same with Figure [51 All other decoupled systems have the 
same bifurcation patterns but without the pair of persistent eigenvalues. For the 
entire spectrum of the 3D linear NS operator, the limiting picture is the same with 
Figure HH as e ^ 0+; and the spectrum is the same with Figure [T2l when setting 
e = |28j . It seems that the Unstable Disk Theorem |10j of the 2D linear Euler 
case is still valid: \k + jp\ < \p\ for some j, implies that there is an eigenvalue of 
positive real part; while \k + jp\ > \p\ for all j, implies that there is no eigenvalue 
of positive real part. 



7.2. The ABC Fixed Point. In this case, the periodic domain is the cube, 
i.e. a = /3 = 1. The ABC flow is given specifically by 

fll — A sin mx^+C cos mx2, — B smmxi+Acosmxs^ fij = Csinmx2+-Bcosmxi, 

where m is a positive integer, and (A, B, C) are real parameters. In terms of Fourier 
modes: Let p — (m,0, 0), q = (0,to,0), and r = (0,0, m), then the ABC flow is 
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Figure 26. The same as Figure [24l except e = 10 ^, a = 1208, 
and T = 2 X IO-^tt. 



given by 



2 



1 



c.*(g) = -C, ^(-g) = --C, 



2i 



The spectral equation for the hnear 3D NS operator at the ABC flow is then given 

by 



\uji{k) = -e\k\^ujt{k) - ^ 

k=k+k 



ran 

(fc)] 



+ |fc| fc,„cj„(fc)[e£ 

mn m] 
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Figure 27. The same as Figure [24l except e = 10"^, a = 1208.2, 
and T = 2 X IO^tt. 



Calculating the eigenvalues of the Galerkin truncations of this system becomes 
challenging. Beyond the size {|A:„| < 6, n = 1, 2, 3}, the computing time is too long. 
Below we present some pictures for the Galerkin truncation {|fc„| < 4, n — 1,2, 3}. 
We choose m = I, A = 1.2, B = 0.7 and C = 0.9. When e = 20000, all the 
eigenvalues are negative as shown in Figure [HIl As e is decreased, eigenvalues start 
to jump off the real axis and form vertical lines as shown in Figure [32] when e = 10, 
in contrast to the parabolas in the cases of cat's eye and 3D shear. When e is 
decreased to e — 0.1, many eigenvalues move to the right half plane, i.e. there are 
many unstable eigenvalues as shown in Figure [331 Notice that for the full linear NS 
operator, there should be an infinite tail of negative eigenvalues to the left. When 
e = 0, the eigenvalues of the Galerkin truncation of linear NS are symmetric with 
respect to the real and imaginary axes as shown in Figure [34l When e = 0, the full 
linear Euler operator at the ABC flow has a continuous spectrum similar to that at 
the cat's eye [28] . That is, the continuous spectrum of the linear Euler at the ABC 
flow, in any Sobolev space 77" (T^) where s is a non-negative integer, is a vertical 
band of width 2s(t symmetric with respect to the imaginary axis {A : |Re(A) | < sa} 
as shown in Figure [TH where cr > is the largest Liapunov exponent of the vector 
fleld given by ABC flow. Thus the width of the vertical band is proportional to 
the scale s of the Sobolev space _ff''(T^) [28j . The union of all such bands for all 
integers s > is the entire complex plane. The eigenfunctions of the linear NS 
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Figure 28. The eigenvalues of the hne k = (a, 0,0) are all nega- 
tive when e = 2.0. 



at the ABC flow when e > belong to i/*(T^) for all integers s > 0. As e is 
decreased, the eigenvalues move into the right half plane. The e — > 0"*" limiting 
picture of the eigenvalues of the linear NS at the ABC flow is that the eigenvalues 
are dense on the entire plane, in contrast to the left half plane in the case of cat's 
eye as shown in Figure [THl That is, all the eigenvalues of the linear NS at the ABC 
flow condense into the entire plane - "condensation" . Thus the possible instability 
hinted by the right half band of the continuous spectrum of linear Euler in iJ*'(T^) 
may be realized by real viscous fluids in this case in contrast to the cat's eye case. 

8. Numerical Verification of the Heteroclinics Conjecture for 3D Euler 

Equations 

For 3D Euler equations, one can also pose the heteroclinics conjecture. 

• The Heteroclinics Conjecture: There is a heteroclinic orbit of the 3D Euler 
equations that connects fl (|7.4p and —fl. 
As discussed in the Introduction, the rationality of this conjecture comes from the 
fact that 3D Euler equations have a Lax pair [23]. Below we will verify this con- 
jecture for the Galerkin truncation: |k„| < 1 (n = 1, 2, 3) where k = {aKi, (3k2, k^). 
Even though this is the smallest Galerkin truncation, the dimension of the resulting 
system is still very large. For this Galerkin truncation, the fixed point (|7.4p is still a 
fixed point. The linearized Galerkin truncation operator at this fixed point can be 
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Figure 29. A pair of eigenvalues of the line k = (a, 0,0) jumps 
off the real axis when e = 1.8. 

obtained by the corresponding Galerkin truncation the 3D linear Euler operator. In 
this case, the line segment labeled by fc = (a, 0, 0) (e = 0) has a positive eigenvalue 
A = 0.5792, and the corresponding eigenvector v is given by: 

0, -1) 0.328919 - i 0.246347, ^2(1, 0, -1) = -0.246347 - i 0.328919, 
cj3(l, 0, -1) = 0.230243 - i 0.172443, wi(l, 0, 0) = 0, 
cj2(l, 0, 0) = -0.19583 - i 0.26147, 

cj3(l, 0, 0) = 0.183029 - i 0.137081, cji(l, 0, 1) = 0.328919 - i 0.246347, 
W2(l, 0, 1) = 0.246347 + i 0.328919, ^3(1, 0, 1) = -0.230243 + i 0.172443, 

and all other a;^(fc)'s are zero. Starting from the initial condition 

(8.1) w = cj* + lO^^w , 

where w* is the Fourier transform of the fixed point (|7.4p . the approximate hete- 
roclinic orbit reaches order ~ 10"'^ neighborhood of — w* during the time interval 
[0,29.33]. This approximate heteroclinic orbit is the lower branch of the approxi- 
mate heteroclinic cycle shown in Figure 1351 Notice that the approximate hetero- 
clinic orbit here does not have the extra loop as in Figure [231 When more modes 
are included in the Galerkin truncation, extra loops may be generated. 
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Figure 30. A bubble of eigenvalues of the line k — (a, 0,0) has 
developed when e — 0.0007. 



9. Melnikov Integral and Numerical Simulation of Chaos in 3D 
Navier-Stokes Equation 

Setting e = in the 3D NS (|7.ip , one gets the corresponding 3D Euler equation 
for which one has the following constants of motion: 



E 



\urdx , H : 



T3 



u ■ ndx 



T3 



where E is the kinetic energy and H is the helicity. We will use the constant of 
motion 



G^E-H^ I \u 



- • fldx 



T3 



to build a Melnikov integral for the corresponding 3D Navier-Stokes equation (|7.ip . 
We will try to make use of the Melnikov integral as a measure of chaos and to 
conduct a control of chaos, around the 3D shear flow (|7.4p . The gradient of G in u 
or Q is given by 
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Figure 31. The eigenvalues of the (Galerkin truncation of) linear 
NS at the ABC flow when e = 20000. 

where curl = Vx . The gradient is zero at the 3D shear flow (j7.4|) . We define the 
Melnikov function for the 3D NS (fTTjl as 

M = r I ^^G[^^ + f{t,x)+b5{x)]dxdt 

= r f 2 curr^ (u -n)\ An + fit, x)+bS(x)]dxdt . 

Next we conduct numerical simulations on the Gcilerkin truncation \nn I < 1 
(n — 1, 2, 3). When e = 0, the Liapunov exponent cr = for all the numerical tests 
that we run. This indicates that there is no chaos when e = 0. Often the smallest 
Galerkin truncation is an integrable system |12j [22j . In such a circumstance, the 
Melnikov integral represents the leading order term of the distance between the 
broken heteroclinic orbit and the center-stable manifold of the fixed point. But the 
dimension of the center-stable manifold is large. The zero of the Melnikov inte- 
gral implies that the unstable manifold in which the broken heteroclinic orbit lives, 
intersects with the center-stable manifold. Therefore, there is a new heteroclinic 
orbit which lives in the intersection. Such a heteroclinic orbit does not immediately 
imply the existence of chaos, even though it may lead to some transient chaos char- 
acterized by finite time positive Liapunov exponent (infinite time positive Liapunov 
exponent is zero). To compute the Melnikov integral, we choose the external force 
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Figure 32. The eigenvalues of the (Galerkin truncation of) hnear 
NS at the ABC flow when e = 10. 



and control as follows 



/i = asint cos(xi + ax2) , /2 = /s = , 
5i(a;)=^e^'=- , 4 = ^ = 0, 

K 



where the sum is over the Galerkin truncation. Then the Melnikov integral has the 
expression 



(9.1) 



M = Mi + aJM^ + M^ sm{to + 0) + 6M4 , 
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Figure 33. The eigenvalues of the (Galerkin truncation of) linear 
NS at the ABC flow when e 0.1. 



where 



Ms 
Mi 



. Ms M2 

sm = — = , cos = 



^/Mi + Mi' VMi + Ml ' 

^ Re I ieemnkm{'i\k\~^£nsrksi^r - ^n)i^e{k) > dt , 

k 

M2 = COStRe{i\k\^'^ei„inkm{i\k\^'^ensrksUJr- i^n)},,,^Qj^)dt , 



(a, 0,1) 

<y —00 

/ + 00 
sint Re {i\k\^'^eimnkm{i\k\^'^ensrksUJr - t^n)}j.^(„_o,i) ' 

/+00 
{i\k\ '^SnsrksLOr - UJn)} dt , 
-00 I, 



where the sum is over the Galerkin truncation, all the integrals are evaluated along 
the lower heteroclinic orbit in Figure [35] for the time interval [—29.33/2,29.33/2], 
rather than (— cx),oo). Direct numerical computation gives that 

Ml = 645.7 , M2 = 1.581 , M3 = , M4 = 47.86 . 
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Figure 34. The eigenvalues of the (Galerkin truncation of) linear 
NS at the ABC flow when e = 0. 



When 6 = (no control) , M has roots when 

\a\ > 408.4 . 

When e = 10"^ 6 = 0, and T = 4 x IO^'tt, we find that 



300 

: 0.8 X 10- 



a = 400 

(T = 1.3 X 10" 



420 

: 0.8 X IQ- 



600 

: 4.9 X 10- 



800 

: 4.9 X 10-4 



Around a — 400, a has a jump of one order which seems to be in agreement with 
the Melnikov prediction. However, when a — 100, a — 3.6 x 10^'' which shows that 
the Melnikov prediction is not very effective. We do not know the specific reason. 

10. Melnikov Integral and Control of Chaos in 3D Navier-Stokes 

Equation 

Now we turn on the control {b ^ 0). When 

b = -Ml/Mi ~ -13.5 , 

the Melnikov integral M (|9.ip has roots for any a. 

When e = 10-^ b = -13.5, and T = 4 x IO^^tt, we find that 

a 1 a = 10 a = 100 a = 200 

cr = 0.098 a = 0.125 a = 0.095 a = 0.083 . 
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Figure 35. An approximate heteroclinic cycle of the Galerkin 
truncation: < 1 (n = 1,2,3) of the 3D Euler equations. 

Thus under the control, chaos exists even when a = 1. Also changing the value of 
b, a does not change dramatically. Therefore, our control prediction above may not 
be very effective. 

11. Numerical Verification of the Heteroclinics Conjecture for a Line 

Model 

Returning to the 2D Navier-Stokes equation (|2.1|) . numerical simulations on 
large Galerkin truncations are still challenging to the current computer ability. Here 
we will study a simple line model 12 obtained by a special Galerkin truncation [12j . 
Let p — (1, 0) and k = (0, a), the line model is given by the Galerkin truncation: 

{±p, ±{k + np), Vn e Z} 

We will work in the invariant subspace where w^'s are real- valued. The governing 
equation of the line model is 

(11.1) +e[-(n2 + a^)oj„ + F„ + feA„] , 

(11.2) = - ^ A„_i^„tj„_iw„ + e[-cj* + + , 
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where LUn — ^k+np^ ^* — similarly for F and A as the Fourier transforms of / 
and (5, and 

1 



An = '2A{p, k + np) — a 



1 



'n? + Oi^ 

An-i.n = 2A{k + (ri — l)p, k + np) = a 



1 



{n — 1)2 + 11^ + 



For the hne model, verification of the heteroclinics conjecture is relatively easier. 
First of all, for the line model (e = 0), it can be proved that the fixed point 
= 2cosa:i has a 1-dimensional local unstable manifold . The basic idea of the 
proof is that one can apply the Riesz projections to the spectrum of the linearized 
line model operator at the fixed point, and the nonlinear terms have bounded 
coefficients so that they are quadratic in a Banach algebra. For the full 2D Euler 
equation, the difficulty lies at the fact that the nonlinear term is not quadratic in 
a Banach algebra. 

Denote by E the 1 co-dimensional hyperplane 

S I w(i^o) = 0} . 

We have the corollary of Theorem 13.11 

Corollary 11.1. Assume that W" H S 7^ 0; then the heteroclinics conjecture 
is true, i.e. there is a heteroclinic orbit that connects f2 = 2cosa;i and — fJ. 

For any truncation (|n| < N) of the line model, we first calculate the unstable 
eigenvector. Then we track the heteroclinic orbit with the initial condition provided 
by the unstable eigenvector. Numerically exact heteroclinic orbit is obtained for 
any N {\n\ < N). That is, for any N {\n\ < N), it can be verified numerically that 

For \n\ < 32, the heteroclinic orbit is shown in Figure [36l In comparison with the 
full 2D Euler equation, the hyperplane E here is only 1 co-dimensional. This is 
the simplest nontrivial case to study the intersection n E. We also conduct 
calculations on the Liapunov exponents. When e = 0, \n\ < 32, a — for all the 
computed time intervals, which means that there is no chaos. This is true for any 
N [n < N) and any computational time interval. This indicates that the line model 
may be integrable when e = 0. From these facts, it is clear that the line model is a 
good starting point for a rigorous analysis. For instance, it is hopeful to make the 
Melnikov integral theory rigorous. 



12. Melnikov Integral and Numerical Simulation of Chaos in the Line 

Model 

For the line model, the kinetic energy and enstrophy are still invariants when 
e = 0. Choosing the same external force (|5.ip and control (|6.ip . we have the 
Melnikov integral which is the same with that of 2D NS except that the Fourier 
modes summation is over the line model. 



(12.1) 



Mo = Ml + aJMj + Ml sin(io + 9) + bMc 
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Figure 36. Numerically exact heteroclinic orbit of the line model 
(e = 0) for \n\ < 32. 



For the truncation < 32, we evaluate these integrals along the heteroclinic orbit 
in Figure 1361 and obtain that 

Ml = -6.0705, M2 = -0.10665, M3 = 0, = 11.9728 

For the case of no control (6 = 0), when 

\a\ > 56.92 

the Melnikov integral M has roots. 

We conduct some numerical simulations on the (transient) chaos. When e — 
10-^, and T = 2 X IO^tt, we find that 

a < 200 a = 400 

cr < cr = 7.2 X 10"^ . 

According to the roots of the Melnikov integral, when \a\ > 56.92, the broken 
heteroclinic orbit may re-intersect with certain center-stable manifold. According 
the above Liapunov exponent result, transient chaos is generated when a = 400 
which may be due to the generation of new heteroclinic cycles. As can be expected, 
the Melnikov prediction performs better here for the line model. When e = 0, there 
is no chaos. When e > 0, we see that transient chaos appears when \a\ is greater 
than certain threshold which is in the range \a\ > 56.92 predicted by the Melnikov 
integral. 
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13. Melnikov Integral and Control of Chaos in the Line Model 

Now we turn on the control (6 7^ 0) . When 

b = -Mi/Mc = 0.50702425 

the Melnikov integral M has roots for any a. To test the effectiveness of the 
control, we set b to the above value and conduct some numerical simulations on the 
(transient) chaos. When e = 10"^, b = 0.50702425, and T = 2 x IO^tt, we find that 

a = 1 a 10 a = 50 a 200 a = 400 

CT = -12.6 X 10"-* (7 = 2x10-4 (7 = 0.2x10-4 cr = 2 x 10^^ cr = 0.87 x 10"^ . 

The control clearly enhanced the transient chaos. The control effectively pushed 
the threshold of a backward from 400 to 10 for the generation of transient chaos. 
This shows that the Melnikov integral control performs better here for the line 
model. 

14. Numerical Verification of the Heteroclinics Conjecture for a Two 

Lines Model 

To gain an understanding of the effect of the other modes np {\n\ > 2) on the 
line model, we introduce the two lines model which is the Galerkin truncation: 

{{kuk2), |fc2|<l}. 

We also work in the invariant subspace where cok's are real- valued. 

For the two lines model, one can derive the governing equations in the physical 
variables. Let 

n = uj{t, x) + e'"yq{t, x) + e-'^yqit, x) , 

where to is real-valued, q is complex-valued (the Fourier transform of q is real- 
valued), and 

/■27r 

/ uj(t, x)dx = . 

Let 

f + bS = r]it, x) + e'"^F(t, x) + e-'^^Fit, x) , 

where 77 is real- valued, F is complex- valued (the Fourier transform of F is real- 
valued), and 

p27T 

/ ri{t, x)dx = . 
Jo 

Substituting the above expressions into the 2D NS (|2.ip . and ignoring the terms 
involving e*^"^ and e"*^"**, one gets the two lines model in the physical variables, 

(14.1) idtq + a[{d^uj){dl-a')-^ -{d-^u)\q = ie[{dl-a')q + F] , 

(14.2) dtio + iad, [q{dl - a^r^q - q{dl - a^Y^q] = e [dluj + r]\ . 
Introducing 6 — d^^us, ip — {d^ — a^)^^q, and h = d^^rj, one gets 

(14.3) idtq + aiipdle-9q) = ie[{dl-a^)q + F] , 

(14.4) dte + ta{q'f-qip)^e[dle + h], 

(14.5) {dl-a')ip = q. 
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When e = 0, then Kinetic energy and enstrophy 

Eq^ [9^ + 2a^\ipf + 2\da,.f\^]dx , Ei= [lo^ + 2\q\^]dx , 
Jo Jo 

are still invariants. 

Denote by E the hyperplane 

Y, = {lu \ LOk = , whenever k2 = 0} ■ 

For the two lines model, when e = 0, existence of a local unstable manifold for the 
fixed point = 2cosa;i is an open problem due to the fact that the coefficients of 
the nonlinear terms are not bounded, i.e. the nonlinear terms are not quadratic in 
a Banach algebra. Then we have the corollary of Theorem 13.11 

Corollary 14.1. Assume that the fixed point fl = 2 cos xi has a 1- dimensional 
local unstable manifold , and W" n E ^ 0; then the heteroclinics conjecture is 
true, i.e. there is a heteroclinic orbit that connects CI = 2cosa;i and — il. 

For any truncation < N) of the two lines model, we first calculate the 
unstable eigenvector. Then we track the heteroclinic orbit with the initial condition 
provided by the unstable eigenvector. For \ki\ < 2, numerically exact heteroclinic 
orbit is obtained. That is, it can verified numerically that 

ty" n E 7^ . 

Figure [571 shows the numerically exact heteroclinic orbit. For |fci| < 4, 

Distance {W'',T.) w 0.0086 . 

Figure [38] shows the approximate heteroclinic orbit. For \ki\ < 16, 

Distance (ty",S) « 0.012 . 

Figure [39] shows the corresponding approximate heteroclinic orbit. Unlike the line 
model, here we do not always get numerically exact heteroclinic orbits. This is due 
to the influence of the modes {ki, 0). 

15. Melnikov Integral and Numerical Simulation of Chaos in the Two 

Lines Model 

When e = 0, there is very weak chaos for the computational interval t E 
[0,4 X IO-^tt]: 

|n|<4 |n|<8 |n| < 16 

cr = 9.5x10-4 (7 = 8.2x10-4 cr = 8.5x10-4. 

For the two lines model, the kinetic energy and enstrophy are still invariants 
when e = 0. Choosing the same external force ()5.1|) and control (|6.ip . we have the 
Melnikov integral which is the same with that of 2D NS except that the Fourier 
modes summation is over the two lines model, 

(15.1) Mo = Ml + a^jMi + Ml sin(io + 0) + bMc . 

For the truncation |fci| < 16, we evaluate these integrals along the heteroclinic orbit 
in Figure [39l and obtain that 

Ml = -4.9 , M2 = -0.0948 , M3 = , Mc = 12.2498 . 
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Figure 37. Numerically exact heteroclinic orbit of the two lines 
model for |fci| < 2. 



For the case of no control {b = 0), when 

\a\ > 51.688 , 

the Melnikov integral M has roots. We conduct some numerical simulations on the 
(transient) chaos. When e = 10^^, and T = 2 x IO^tt, we find that 

a < 200 a = 400 

cr < (7 = 1.27 X 10-2 . 

According to the roots of the Melnikov integral, when \a\ > 51.688, the broken 
heteroclinic orbit may re-intersect with certain center-stable manifold. According 
the above Liapunov exponent result, strong transient chaos is generated when a = 
400 which may be due to the generation of new heteroclinic cycles. Thus the result 
is almost the same as that of the line model. As can be expected, the Melnikov 
prediction also performs well here for the two lines model. When e = 0, there is 
weak transient chaos. When e > 0, we see that strong transient chaos appears when 
|a| is greater than certain threshold which is in the range \a\ > 51.688 predicted by 
the Melnikov integral. In comparison with the line model, here the transient chaos 
seems very strong. On the other hand, when e = 0, the two lines model here still 
has weak chaos, in contrast to the fact of no chaos at all for the line model. 



48 



YUEHENG LAN AND Y. CHARLES LI 








Figure 38. Approximate hetero clinic orbit of the two lines model 
for |fci| < 4. 

16. Melnikov Integral and Control of Chaos in the Two Lines Model 

Now we turn on the control (6 7^ 0). When 

b = -Mi/Mc = 0.4 , 

the Melnikov integral M has roots for any a. To test the effectiveness of the 
control, wc set b to the above value and conduct some numerical simulations on the 
(transient) chaos. When e — f 0^^, b = 0.4, and T — 2 x IO'^tt, we find that 

a < 400 a = 500 a = 600 a = 700 a = 800 

C7 < (7 = 0.97 X 10-3 cr = 3.0 X 10"^ ct = 4.6 x lO'^ cr = 5.0 x IQ-^ . 

The control seems to tame the transient chaos, in contrast to the line model. The 
control pushed the threshold of a forward from 400 to 700 for the generation of 
strong transient chaos. 

17. Conclusion and Discussion 

Through a combination of analytical and numerical studies, we now have a 
better understanding on the zero viscosity limit of the spectra of linear NS oper- 
ators. We proposed and numerically studied the so-called hcteroclinics conjecture 
for both 2D and 3D Euler equations. We also proposed the Melnikov integral as a 
tool for predicting and controling chaos. 



NAVIER-STOKES AND EULER EQUATIONS 



49 




Figure 39. Approximate hetero clinic orbit of the two lines model 
for |fci| < 16. 



Our numerical verification on the heteroclinics conjecture was limited by reson- 
able computing time. We realized that increasing the size of the Galerkin trunca- 
tions can quickly reach the limit of our computer ability. We did not try to utilize 
today's supercomputer due to the fact that Galerkin truncations arc essentially 
singular perturbations of Euler equations. We believe that analysis is the key to a 
better understanding of the heteroclinics conjecture. 

We realized through our numerical simulations that Liapunov exponent per- 
forms very well as a measure of (even transient) chaos. Our Melnikov prediction 
for NS equations is of course not as rigorous and effective as that for sine-Gordon 
system. In fact, it is a rough indicator for predicting and controling chaos. Nev- 
ertheless, we believe that the Melnikov integral theory for NS equations has a lot 
of potential especially in the current circumstance that there is no effective tools 
dealing with chaos in NS equations. 

We also believe that both the line and the two lines models have great potential 
in future analytical studies on modeling the dynamics of 2D NS equations. 
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(18.1 



18. Appendix A: Melnikov Integral and Control of Chaos in a 

Sine-Gordon Equation 

Consider the sine-Gordon equation |15| 
9 



Utt 



16 



+ sm u 



auf + [l + bS{x)] cos t sin u 



which is subject to periodic boundary condition and odd constraint 

(18.2) u{t,x + 2tt) = u{t,x) , u{t,-x) ^ ~u{t,x) , 

where u is a real- valued function of two real variables (t, x), e is a small perturbation 
parameter, a > is the damping coefficient, bS{x) is the spatially localized control, 
S{x) is an even and 27r-periodic function of x, and b is the control parameter. The 
system (jl8.1l) is invariant under the transform u — )■ —u. 

The natural phase space for (|18.ip is {u,ut) e x i?" (n > 0) where i?" 

is the Sobolev space on [0,27r]. Let P be the Poincare period-27r map of (|18.1|) in 
fjn+i ^ jjn •\Yithout the control (6 = 0), we have the following chaos theorem |15j 

Theorem 18.1 ( [15j ). There exists a constant gq > 0, when e is sufficiently 
small, for any a € [i^QOi Qq] there exists a symmetric pair of homoclinic orbits h± 
(7i„ — ~h^) asymptotic to (u,Ut) = (0,0). In the neighborhood of h±, there exists 
chaos to the sine- Gordon equation U8.1\) in the following sense: There is a Cantor 
set S of points in H"~^^ x (n > 0), which is invariant under an iterated Poincare 
period-2'K P^ for some K. The action of P^ on S is topologically conjugate to the 
Bernoulli shift on two symbols and 1. 

In the product topology, the Bernoulli shift has the property of sensitive de- 
pendence upon initial data - the signature of chaos. 

When we turn on the control (5^0), we hope to find values of b such that the 

annihilated or less chaotic, enhanced 



chaos in Theorem II 8. II is controlled (tamed 
- more chaotic). Our main tool is the Melnikov function. To build such a function, 
we need results from integrable theory. When e = 0, the fixed point (u = 0) of the 
sine- Gordon equation (|18.ip has a figure eight connecting to it 



(18.3) 



±4 arctan 



V7 



sech r sin x 



where r — ^{t — to) and to is a real parameter. Along this figure eight, a Melnikov 
vector has the expression |15| : 



dFi 7n 

(18.4) — — = ± = sech rtanhrsinx 

dut 12V2 



9 7 ,2 ■ 2 ' 
— -I- — sech r sm x 
16 16 



where Fi is a constant of motion. When e 7^ 0, the Melnikov function for p8.ip is 
given by [15) : 

n+00 r2n gp^ 



M{to,a,b) = 



dut 



-aut 



(1 + b8{x) 



cos t sin u 



dxdt , 



where u and are given in (jl8.3p and (118. 4p . Using the odd and even property 
of p8.3p and (|18.4p in t and x, we obtain 



(18.5) 



M {to, a, 6) = -aMa + sin to (Mo + bMb) , 
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where 



+00 ^27r Qp 

Ma = I / -7^ — ut dxdt , 
00 Jo out 
f+00 i'2-iT gp_^ ^ 
Mq = — / / — — sin — -=T sin^ u dxdt , 



00 



dut 



/+00 i<2tt qp _ ^ 
/ — — S(x) sin —=T sin^ m dxdt . 
-00 Jo 9ut ' ' V7 

In the phase space x 7?" («■ > 0), {u, Ut) — (0, 0) is a saddle point under the 

Poincare period-27r map of (|18.ip with one-dimensional unstable manifold and 
one-codimensional stable manifold . The Melnikov function eM{to, a,b) is the 
leading order term of the distance between and W. For the entire rigorous 
theory, see [16]. When \aMa\ < \Mo + bMb\, the roots of M are given by 

(18.6) smto 



Near these roots, and W intersect. This leads to the existence of a symmetric 
pair of homoclinic orbits and chaos in Theorem 118. II When |aAfa| > \Mo + bMt,\, 
i.e. 

(18.7) - a\Ma\ - Mo < bMb < a\Ma\ - Mq , 

the Melnikov function is not zero, and we have the following theorem. 

Theorem 18.2. When the control parameter b satisfies 7[ j, the chaos in 
Theorem \18.1\ disappears. 

Proof. When the control parameter b satisfies ()18.7|) . the Melnikov function 
is not zero for any to, and and W do not intersect. Thus the pair of homoclinic 
orbits and the corresponding chaos in Theorem 1 1 8 . 1 1 disappear . □ 



Theorem ll8.2l onlv claims that the chaos in Theorem ll8.1l disappears. This does 
not mean that there is no chaos in the entire phase space x {n> 0). An 

important point here is that by manipulating the localized control bS{x), one can 
change the Melnikov function which leads to the disappearance of the non-localized 
(in x) chaos. The control condition (|18.7p is also interesting: It is not true that the 
larger the control parameter b is, the better the taming is. In fact, when b is large 
enough, the chaos will reappear. 

When \aMa\ < |A/o + bMb\, the Melnikov function P^3)) has roots ^SMi and 
Theorem llS.ll holds. As a function of to, the Melnikov function M has the maximal 
absolute value (the L°° norm), 

M4a,b) ^ a\Ma\ + \Mo + bMb\ , 

for to € [0, 2tt]. M* is the leading order term of the maximal distance between 
and W^. Notice that and intersects near the to given by (|18.6p . So the 
larger the M■^, is, the more violent the chaos is. Thus M^,{a,b) serves as a measure 
of the strength of the chaos. By changing the control parameter b, we can adjust 
the strength M* of the chaos - enhancing or decreasing. 
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19. Appendix B: The Lagrange Flow Induced by a Solution to the 2D 
Euler Equation Is Always Integrable 

It is well-known that 2D Euler equation is globally well-posed [5] [6]. For any 
solution to the 2D Euler equation, let ^' = 'I'(t, xi, X2) be the corresponding stream 
function. Then the Lagrange flow induced by the solution is given by 

(19 1) dx2 _ 

dt dx2 ' dt dxi 

Theorem 19.1. The Lagrange flow \19.1\] induced by a solution to the 2D Euler 
equation is always integrable. 

Proof. Assume that 'i/{t,xi,X2) is not a steady state, i.e. it depends upon 
t (in this case A^f is functionally independent of 5", otherwise, ^(t,xi,X2) would 
be a steady state). Introducing the new Hamiltonian H = '^(9,xi,X2) — 'ip, and 
converting (jl9.ip into an autonomous system 

dxi _ dH dx2 _dH dO _ dH _ dH 

^ ~dr~~'dx^' ~dr ~ dx^ ' d< ^ " ' 'dt~~d9' 

Notice that the vorticity fl — A'^ is another constant of motion of (|19.2p besides 
H: 

^A^ = deA^ - d,,A^d,,^ + 9,, A*a,,* = . 
dt " " 

Since A^f is independent of ip, A^f and H are functionally independent. Thus 
(119. 2p is integrable in the Liouville sense. In the case that ^' is independent of t 
(i.e. a steady state), then (|19.1[) is an autonomous system, thus also integrable in 
the Liouville sense. □ 

A common way to obtain steady states of 2D Euler equation is by solving 

A*-/(*) 

where /(^') is an arbitrary function of 5*, i.e. A^" and 5* are functionally dependent. 
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